Shell Model Monte Carlo method in the pri-formalism 
and applications to the Zr and Mo isotopes 



in 

O 

o 

< 

in 



> 
o 

oo 
O 
in 
o 



o 

=5 



X 



C. Ozen 1 and D.J. Dean 2 
1 Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996 and 
2 Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831 

(Dated: February 5, 2008) 

We report on the development of a new shell-model Monte Carlo algorithm which uses the proton- 
neutron formalism. Shell model Monte Carlo methods, within the isospin formulation, have been 
successfully used in large-scale shell-model calculations. Motivation for this work is to extend the 
feasibility of these methods to shell-model studies involving non-identical proton and neutron va- 
lence spaces. We show the viability of the new approach with some test results. Finally, we use 
a realistic nucleon-nucleon interaction in the model space described by (lpi/2, Ogg/2) proton and 
(lds/2i 2si/2, 1^3/2,0^7/2,0/111/2) neutron orbitals above the 88 Sr core to calculate ground-state en- 
ergies, binding energies, B(E2) strengths, and to study pairing properties of the even-even 90_104 Zr 
and 9 ~ 



3 Mo isotope chains. 



I. INTRODUCTION AND MOTIVATION 

The Shell model Monte Carlo (SMMC) method 
was developed as an alternative to direct diagonaliza- 
tion in order to study low-energy nuclear properties. It 
was successfully applied to nuclear problems where large 
model spaces made diagonalization impractical. In the 
canonical SMMC approach, one calculates the thermal 
expectation values of observables of few-body operators 
by representing the imaginary-time many-body evolution 
operator as a superposition of one-body propagators in 
fluctuating auxiliary fields. Thus, one recasts the Hamil- 
tonian diagonalization problem as a stochastic integra- 
tion problem. 

In this paper, we report on the development of an 
SMMC approach in the pn-formalism. This implementa- 
tion of SMMC enables one to treat shell-model Hamilto- 
nians that are not isospin invariant in the model space, 
or for which different model spaces are used for protons 
and neutrons. In the following, we will use the abbre- 
viated form SMMCpn, to distinguish the approach dis- 
cussed here from the original one. We note that special 
features of the pairing+quadrupole interaction enabled a 
special implementation of SMMC in non-degenerate pro- 
ton and neutron model spaces for calculations in rare 
earth nuclei 0, || . The method presented in this work 
is general and may be used for realistic Hamiltonians, as 
well those of a more schematic variety. 

As a first novel application of the new implementa- 
tion, we perform shell-model calculations for the even- 



even 



90-104 



Zr and 



92-106 



Mo isotopic chains. Initial ex- 



perimental studies |6j indicated that nuclei in this re- 
gion have very large deformations, and that the tran- 
sition from spherical shapes to highly deformed shapes 
occurs abruptly: 96 Zr is rather spherical, while 100 - 104 Zr 
nuclei are well deformed with a quadrupole deformation 
parameter of /?2 = 0.35 0. Furthermore, the spherical- 
to-deformed transition is more abrupt in the Zr isotopes 
than in the nearby elements Mo, Ru, and Pd. Generator- 
coordinate mean-field calculations in this region Q are 



able to reproduce the shape transitions with particular 
Skyrme interactions. Furthermore, the region exhibits 
significant shape-coexistence phenomena [9lll0| . 



The history of shell-model applications in this mass 
region goes back to the 1960s with model spaces built 
on 88 Sr or 90 Zr cores [H IH E H. Gloeckner 
used an effective interaction built on a 88 Sr core with a 
model space consisting of the orbitals ir : (lpi/2, 0g 9 / 2 ), 
v : {X d^ip . 2s 1 / o ). Other studies used larger model 
spaces [l(| E3, 13 with varying effective interactions and 
truncation schemes. Holt et. al. 0] derived a realis- 
tic effective interaction using many-body perturbation 
techniques |2(j in the model space n : (lpi/2, 0*79/2), 
v : (lrf 5 /2)2s 1 /2, ld 3 / 2 > 037/2>0/i u / 2 )- This effective in- 
teraction was based on the realistic nucleon-nucleon CD- 
Bonn potential ;2lj , and shell-model diagonalization cal- 
culations were carried out for the low-lying spectra of 
the Zr isotopes with neutron numbers from iV = 52 to 
N = 60. Their results showed reasonable agreement with 
experimental spectra. In this article, we use a slightly 
modified version of this realistic effective interaction to 
explore Zr and Mo nuclei through N = 64. 



In Sec. [HJ we give an outline of the SMMC method 
with an emphasis on the differences in the SMMCpn im- 
plementation when compared to the isospin-conserving 
implementation. Then, in Sec. IIII Al we will demon- 
strate the utility of the new approach by a comparison 
of various numerical results obtained using the SMMCpn 
technique for a few /p-shell nuclei to those calculated by 
direct diagonalization and earlier SMMC studies. Calcu- 
lations for the Zr and Mo isotopes, which are presented 
in Sec. IIII 51 were carried out in the same model space as 
in Holt et. al. 0], using a slightly modified interaction 
[22|. We show results for ground-state energies, binding 
energies, B(E2) strengths, and BCS-like pairing corre- 
lations for the Zr and Mo isotope chains. We conclude 
with a perspective on this avenue of research. 
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II. FORMALISM 



where the metric of the functional integral is 



In the SMMC method, we calculate expectation values 
of operators within a thermal ensemble of particles whose 
interactions are governed by the Hamiltonian H of the 
system. (A zero-temperature formalism also exists but 
will not be discussed here.) The canonical expectation 
value of an operator A at a temperature T is given by 



(X) 



Tr[UX] 



(1) 



where the partition function of the system is given by 
Z((3) — TrU, j3 — 1/T is the inverse temperature (with 
units MeV -1 ), and the many-body evolution operator is 
U = e~@ H . The quantum-mechanical trace of an opera- 
tor is defined as 



TrA = 



A 



(2) 



where the sum runs over all many-body states in the 
Hilbert space. For nuclear calculations, the number of 
valence particles is usually limited, so that number pro- 
jection becomes important. The original SMMC method 
preserved isospin within the same neutron and proton 
model spaces. We discuss in the following how to im- 
plement number projection when the isospin quantum 
number is broken. Note that in the limit of (3 — > oo, 
we may evaluate ground-state properties of the nuclear 
systems. 

In the following, we will consider Hamiltonians that 
have at most two-body terms. Any such Hamiltonian 
can be cast into a quadratic form: 



H = ^e a p a + ^^A Q /5 2 , 



(3) 



2tt 



and the Gaussian weight is given by 



Q = e -^A/3E Qi „|A»|^(r„) 



The one-body evolution operator is written as 



(5) 



(6) 



(7) 



where we note the dependence on the auxiliary fields 
cr(r„). This time-ordered product means that this for- 
malism yields a path integral in the fields a. The lin- 
earized one-body Hamiltonian for the time slice r n is 
given by 



(8) 



with s = ±1, if A > 0; or s = ±i, if A < 0. Note that 
because the various p a need not commute, (3.18) is accu- 
rate only through order A/3 and that the representation 
of e~ A/3h must be accurate through order A/3 2 to achieve 
that accuracy. 

The thermal expectation values can be expressed as 
the ratio of path integrals in fluctuating auxiliary fields, 

- Tr[Oe-n /P[q]G CT (6U. 
{ ' Tr[e-e&] JV[a]G^ a W 

where the following definitions are used: 



where s a is the energy of single-particle level a, and 
the operator p a is a a one-body density operator of the 
form a^a. Details are given in Q on how to transform 
one-body operators with quantum numbers {n, I, j, j z , t z } 
(where n is the principal quantum number, I is the orbital 
momentum, j is the total angular momentum, and j z is 
its projection, and t z = ±1 for protons and neutrons) to 
the form shown in Eq. J3J. 

At the heart of the SMMC method lies the linearization 
of the imaginary-time many-body propagator. Since, in 
general, \p a ,pp] ^ 0, we must split the interval /3 into 
N t "time slices" of length A/3 = 0/N t We apply the 
Hubbard- Stratonovich transformation |23ll24l ] to the two- 
body evolution operator at each time slice. In compact 
notation, the partition function can be written as: 



Z(p) = TrU = Tie 



-0H 



Tr 



-A0H 



N t 



D[<T}G(a)TrU a 



(4) 



= TrU a 



(0)a = 



Tr[Oc/ CT 
TrU*? 



(10) 



In order to use Metropolis Monte Carlo sampling meth- 
ods [2!| , we need to define a positive-definite weight func- 
tion, 



W a = G (T \&\, 
so Eqn. El can now be rewritten as 



(o> = 



fV[a]W a (d) a Q a = ((6) Mw 



where 



(11) 



(12) 



(13) 



3 



is the sign of the partition function. 

The description above shows how one may transform 
the shell model into a problem of quadrature integration. 
Objects, £ CT and {X)„, in the integrands are of one- body 
nature and are represented by N s x N s dimensional ma- 
trices where N s is the number of the single-particle levels 
in the valence space. The path integrals in the auxiliary 
fields are evaluated by performing a Metropolis random 
walk in the held space. Thermodynamic expectation val- 
ues are given as the ratio of two multidimensional inte- 
grals over the auxiliary holds. The dimension D of these 
integrals is of order N^N t , which can exceed 10 5 for the 
problems of interest in this paper. 

Note that the Monte Carlo sign problem enters calcu- 
lations when any of the A Q matrix elements is positive. 
Realistic shell-model interactions always have such terms; 
a special case is the pairing-plus-quadrupole Hamiltonian 
which has no sign problems. 

If the proton and neutron valence spaces are identi- 
cal and the Hamiltonian H is isospin-symmetric, then 
the Hamiltonian can be cast into a quadratic form which 
respects this symmetry explicitly. In that case, it is pos- 
sible to form linear combinations of density operators 
that separately conserve the neutron and proton num- 
bers. In the isospin formulation (as done in the origi- 
nal SMMC studies) proton and neutron many-body wave 
functions are represented by Slater determinants, and the 
ensuing one-body propagator factors into two propaga- 
tors as well, one for protons and another for neutrons. 
The canonical traces are then calculated by applying the 
number-projection operator to obtain the desired proton 
and neutron numbers. In contrast, we will employ the 
shell-model Monte Carlo method to Hamiltonians that 
are not necessarily isospin invariant or for which proton 
and neutron valence spaces are different. A relaxation 
of the isospin symmetry in the quadratic forms becomes 
essential in order to employ the SMMC method in such 
cases. In the pn-formalism, proton and neutron valence 
spaces are no longer distinguished from each other; in- 
stead, we consider a single valence space containing both 
proton and neutron orbitals. In this way, the density 
operators in the one-body Hamiltonian h{a) inevitably 
mix protons and neutrons, and as a consequence their 
respective expectation values will fluctuate from sample 
to sample. The canonical trace is then retrieved by em- 
ploying projection operators to fix the total number of 
particles A and the z-component of the total isospin T z . 
This implementation represents the major difference be- 
tween the original SMMC and the SMMCpn techniques. 

Projection operators for hxed A and T z are given by 



Pa 



2tt 



2tt 



and 



277 



(16 



r> I — -iBT z i6T, 



(14a) 



(14b) 



respectively. In the discrete Fourier representation, we 
make the substitution. 



2 * d9 

(i 2-7T 7o 27T 



N 3 N Tz 



N s (N Tz + 1) 



where the quadrature points are given by <f> m = 2irm/N s 
and 8 n — 2irn/(NT z + 1), and Nt z is the number of values 
T z can take. As an example, the canonical trace of the 
one-body propagator U a can be obtained by acting with 
both projection operators on the grand-canonical trace, 
TrU a = det(l + U ff ): 



Tr 4 t, U a 



r_ + 1) Z - 



N S (N, 



1 y 

h, + 1) ^ 



N s (N Tz + 1) 
x det(l + e^e ia '^V a ) , 



e -i<pA e ~iST z 



e -iM e -i9T z x 



m.n 



(16) 



(17) 



where the boldface symbols are used for the matrix rep- 
resentation of the operators, for example 



(18) 



-1 










A typical difficulty in the SMMC applications is due 
to a sign problem arising from the repulsive part of the 
realistic interactions. In the case of realistic interactions, 
a straightforward application of Eqn. lllfl to obtain a 
positive definite weight will introduce a highly fluctuating 
integral. This will give rise to expectation values with 
very large fluctuations. In order to avoid this situation, 
we adopt a practical solution |26( to the sign problem by 
breaking the two-body interaction into "good" (without 
a sign problem) and "bad" (with a sign problem) parts: 
H = Hq + Hb- Using a parameter g, we then construct 
a new family of Hamiltonians H(g) = Hq + qHb which 
are free of the sign problem for non-positive values of 
g. The SMMC observables are evaluated for a number 
of different g values in the interval — 1 < g < 0, and 
the physical values are thus retrieved by extrapolation to 
g = 1. We use polynomial extrapolations and choose the 
minimum order that gives \ ~ 1. In our calculations, 
most of the extrapolations are linear or quadratic. In the 
extrapolation of (H) , the variational principle imposes a 
vanishing derivative at the physical value g = 1. A cubic 
extrapolation in this case typically gives the best results. 

Another problem encountered in applying the SMMC 
methods concerns efficiency of the Metropolis algorithm 
in generating uncorrelated field configurations. Rather 
than sample continuous helds, where decorrelated sam- 
ples are obtained after only very many (of order 200) 
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Metropolis steps, we approximate the continuous inte- 
gral over each a a (r n ) by a discrete sum derived from a 
Gaussian quadrature In particular, the relation 

daf(a)e Al3x ^ (19) 

-00 

is satisfied through terms in (A/3) 2 if 

f(a) = I [8{a - a ) + 5(a + a ) + 46(a)] , (20) 
6 

where <t = (3/AA/3) 1 / 2 . In the SMMCpn algorithm, 
we find that samples are well decorrelated after only a 
few (typically less than ten) Metropolis steps using these 
descretized fields. 

We describe in the following section our initial results 
using the SMMCpn method. 

III. RESULTS 

A. Comparisons with direct diagonalization and 
previous SMMC calculations 

In this section, we present a number of test cases that 
validate the SMMCpn approach. For this purpose, we 
first carried out calculations on a few sd-shell nuclei using 
a quadrupole-plus-pairing interaction which is free of the 
sign-problem. This interaction can be written as 

V = - X Q ■ Q - gP {0 ^ ■ P {0A) (21) 

where 

q = iEtf-H^^M !. ® ^] J=1 ' T=0 (22) 

v ab 

and 

a 

The strength of the interaction terms were chosen as 
X = 0.0260 MeV^fm 2 and g = 0.212 MeV. We adopted 
the standard USD [27j single-particle energies. The term 
Vws m the equation above is the central part of a Woods- 
Saxon potential with parameters given in [2c|. The 
SMMCpn calculations were performed at /3 = 2 MeV -1 
with N t = 128 time slices (A/3 = 1/64 MeV" 1 ), and 
each calculation involved 2500-3000 uncorrelated sam- 
ples. Note that typical isospin-conserving SMMC cal- 
culations only require A/3 = 1/32 MeV -1 for similar 
convergence. A comparison of our results with the di- 
rect diagonalization values (obtained by running the AN- 
TOINE |29| code) is given in Table D Results are com- 
patible within the internal heating energy and statistical 
errors of the thermal SMMC calculations. 



TABLE I: Comparison of the ground-state energies (in MeV), 
as calculated by SMMCpn and the ANTOINE results. A 
quadrupole-plus-pairing interaction is used. 



Nucleus E (ANTOINE) 


E (SMMCpn) 


24 Mg 


-39.28 


-38.68 ± 0.27 


26 Mg 


-43.58 


-42.75 ± 0.70 


22 Ne 


-30.23 


-29.53 ± 0.41 


28 Si 


-49.25 


-48.80 ± 0.37 



We also tested the viability of the new implementa- 
tion with the utilization of the extrapolation method de- 
scribed above. For this purpose, a few /p-shell nuclei 
were chosen and the modified Kuo-Brown KB3 residual 
interaction [3(| was used. We calculated the ground-state 
energies, total B(M1), B(E2), and the Gamow- Teller 
strengths of the nuclei and compared our results to those 
obtained by exact diagonalization |3l| and those obtained 
by the isospin SMMC The SMMCpn calculations 

were performed at /3 = 2 MeV -1 with A/3 = l/64MeV"\ 
and each calculation at each of six values of the extrapo- 
lation parameter g involved 8000-9000 uncorrelated sam- 
ples. We used a quadratic extrapolation for the total 
B(Ml) and B(GT+) strengths, while for the total B(E2) 
strengths, a linear extrapolation was more reasonable. 
The ground state energies employed cubic polynomials 
subject to the constraint d(H) / dg\ g= \ = due to a vari- 
ational principle that (H) obeys. In all the cases, the 
errors were conservatively adopted from a quadratic ex- 
trapolation. 

In TablellTI ground-state energies and B(E2) strengths 
are tabulated. In all cases, the energies agree strikingly 
well within error bars that are reasonable with the inter- 
nal excitation energy of a few hundred keV due to the 
finite-temperature calculations. The B(E2) strength is 
given by 

B{E2) = ((e p Q p + e n Q n ) 2 ), (24) 

where the quadrupole operator is defined as Q v r n ) — 
J2i r i^2(Si, 4>i). The effective charges were chosen to be 
e p = 1.35 and e n = 0.35, and the oscillator strength is 
given by b = 1.01A 1 / 6 . B(E2) values are also nicely re- 
produced in general, while in the case of Cr, the exact 
result is underestimated by w 25%. 

Table ITTT1 shows a comparison of results for the B(M1) 
and B(GT+) strengths. The B(Ml) strength is defined 
by 

B(M1) = {Q2 m {giU g s s}) 2 ) (25) 

where /xjv is the nuclear magneton. We used the bare 
g factors for angular momentum and spin (gi — 1 and 
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TABLE II: Comparison of exact diagonalization, SMMCpn, 
and isospin SMMC results for valence space energies (in MeV) 
and B(E2) strengths (in e 2 fm 4 ). Typical error bar for energies 
is ±0.6 MeV for SMMCpn and ±0.4 MeV for isospin SMMC 
calculations. 



Nucleus EE E 
exact SMMC(pn) SMMC (iso) 



48 Ti 


-24.6 


-24.4 


-23.9 


48 Cr 


-32.9 


-32.6 


-32.3 


56 Fc 


-66.4 


-66.0 


-65.8 


64 Zn 


-106.3 


-106.5 


-104.8 



Nucleus £B(E2) £B(E2) £ B(E2) 
exact SMMC(pn) SMMC (iso) 



48 Ti 


476 


459 ± 33 


455 ± 25 


48 Cr 


978 


745 ± 40 


945 ± 45 


56 Fc 


1019 


913 ± 55 


990 ± 6 


64 Zn 


1157 


1116 ± 81 


1225 ± 65 



TABLE III: Comparison of exact diagonalization, SMMCpn, 
and isospin SMMC results for B(AI1) (in /in) and Gamow- 
Teller strengths. 



Nucleus 


EB(M1) 


EB(M1) 


EB(M1) 




exact 


SMMC(pn) 


SMMC (iso) 


48 Ti 


10.6 


10.4 ± 4 


10.2 ± 1.2 


48 Cr 


12.0 


12.5 ± 4.4 


13.8 ± 1.7 


56 Fc 


19.4 


22.6 ± 6.4 


20.4 ± 3.0 


64 Zn 


21.6 


22.8 ± 1.2 


23.6 ± 2.2 




Nucleus 


£B(GT+) 


£B(GT+) 


£B(GT+) 




exact 


SMMC(pn) 


SMMC (iso) 


48 Ti 


1.26 


0.89 ± 0.36 


1.13 ± 0.18 


48 Cr 


4.13 


4.35 ± 0.44 


4.37 ± 0.35 


56 Fc 


4.69 


4.02 ± 0.55 


3.99 ± 0.27 


64 Zn 


5.54 


5.66 ± 0.7 


4.13 ± 0.34 



B. Applications to Zr and Mo isotopes 

Calculations for the Zr and Mo isotopes were carried 
out in the same valence space as in Holt et. al. [l9|. 
which is built on the 88 Sr core, but we employed a slightly 
modified interaction which has been previously tested for 
nuclei with small numbers of valence particles in this re- 
gion [23. The model space and single-particle energies 
that were used in the calculations are given in Table ITVl 
r fTg| and references therein). All calculations were per- 
formed at f3 = 2 MeV -1 using N t = 128 time intervals 
with 8500-9500 uncorrelated samples for each extrapola- 
tion parameter g. 

TABLE IV: The model space and single-particle energies used 
in these calculations. 



Protons 




Neutrons 


Orbital Energy (MeV) 


Orbital Energy (MeV) 


0g 9/2 0.90 


0h 11/2 


3.50 


lpi/2 0.00 


0g7/2 


2.63 




ld 3 / 2 


2.23 




2si/ 2 


1.26 




ld 5 /2 


0.00 


e p = 1.8 




e„ = 1.5 


b 


= 2.25 fm 





1. Ground-state energies 

Shown in Fig. ^is the comparison of the expectation 
value of energy (H). Filled circles that are connected 
with dashed lines represent exact diagonalization results 
[33| and in both the Zr and Mo cases, the agreement 
of the SMMCpn values is remarkable. Only for 94 Zr do 
error bars of the SMMCpn result miss the exact diago- 
nalization value slightly. This represents the full test of 
the SMMCpn approach. 



g s = 5.586 for protons and gi = and g s = —3.826 for 
neutrons). The Gamow- Teller strength is defined by 

B(GT+) = (G T G±), (26) 

where the unquenched Gamow- Teller operator is writ- 
ten as G± = Ei^*±. Both B(M1) and B(GT+) re- 
sults agree well with those obtained by direct diagonal- 
ization. Lastly, as in the case of isospin SMMC cal- 
culations, our calculations satisfy the Ikeda sum rule 
B(GT-) - B(GT+) = 3(N - Z) exactly. 



2. Binding energies 

The ground-state energies shown in Fig. correspond 
to the contribution to the nuclear binding energy of the 
interaction of the valence particles among themselves. 

In Fig. 21 we plot calculated and experimental values of 
binding energies with respect to the 88 Sr core. We used 
the following formulae to obtain our binding energies: 

BE( 90+,l Zr) = BE( 90+ "Zr) - BE( 88 Sr) 

-n [BE( 89 Sr) - BE( 88 Sr)] 

-2 [BE( 89 Y) - BE( 88 Sr)] (27) 
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FIG. 1: Ground-state energies of Zr and Mo isotopes. Results 
from direct diagonalization available for lighter isotopes are 
also shown. 



BE( 92+n Mo) = BE( 92+n Mo) - BE( 88 Sr) 
-n [BE( 89 Sr) - BE( 88 Sr)] 
-4 [BE( 89 Y) - BE( 88 Sr)] . (28) 

An inspection of the resulting relative binding energies 
shows that the calculated values (shown by asterisks) de- 
viate from the experimental values (shown by filled cir- 
cles), which display a parabolic behavior. This situation 
is common among calculations using realistic interactions 
derived from AW scattering data [34| . It may be related 
to the absence of real three-body forces in the construc- 
tion of the effective interactions [35| . It is known that a 
given Hamiltonian can always be separated in the form 
H = H m + Hm, where H m is the monopole part, while 
the multipole Hm contains all other terms j3||. Hm 
given by realistic AW interactions takes proper account 
of the configuration mixing; however, the monopole part 
H m fails to produce the correct unperturbed energies. It 
is possible to change the averages of the so-called centroid 
matrix elements to fix this failure without affecting the 
spectrosopy to produce the correct binding energies • 
However, a rigorous treatment of the global monopole 
corrections would deserve a detailed study and thus goes 
beyond the scope of the current work. Instead, to give 
some substance to how a monopole correction may work, 



FIG. 2: Binding energies of the Zr and Mo isotopes. 

we add an overall constant to the diagonal interaction 
elements so that the modified matrix elements are given 
by 

VJ wd (ab, ab) = Vj(ab, ab) + W-^- — '-, (29) 

where n is the number of valence particles. We have 
adopted W — —125 keV to reproduce the binding energy 
of 102 Zr. We plot the effect of this rather naive cor- 
rection in Fig. [3 where modified results, represented by 
diamonds, show much better agreement for both chains 
of isotopes. 

3. B(E2) strengths 

Since the 2\ state is expected to absorb most of the 
total B(E2) strength, the latter can be used as a mea- 
sure of the 0]*~ — 2]*~ spacing, which should reflect a strong 
change with the shape transitions. Shown in Fig. are 
the calculated total B(E2) strengths (open circles) and 
available experimental [38| Of — ^ 2^ values (filled cir- 
cles). Despite the fact that the calculated total strengths 
increase as expected in both isotope chains with the ad- 
dition of neutrons, their numerical values fall somewhat 
below the experimental B(E2; Of — > 2^) values on the 
heavier side of the isotope chains. This failure is quite 
possibly due to our choice of single-particle energy for the 
Q/in/2 which is known to contribute to deformations in 
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FIG. 3: Total B(E2) strengths from the ground state for Zr 
and Mo isotopes. Experimental results are for 0^ — > 2f . 



FIG. 4: Pairing correlations of Zr and Mo isotopes. 



this region 39] ; it may also be due to correlations absent 
due to our choice of model space. We will investigate this 
further in future work. 



4- Pairing correlations 

Pairing correlations among like nucleons is known to 
be important for the ground state properties of the 
even-even nuclei ^(j- These correlations should become 
quenched along the Zr and Mo isotope chains as the tran- 
sition from spherical to well-deformed shapes becomes 
more pronounced. Similar effects were recently investi- 
gated in N = 40 isotones ^lj- We have investigated the 
pairing content of the ground states of the nuclei of in- 
terest using a BCS-like pair operator which is defined for 
neutrons as 



quantity for an uncorrelated Fermi gas is given by 



At = V i/t i/t- 

jm>0 



(30) 



where the sum is over all orbitals with m > 0, and = 
(— l) J+ "Vj_ m is the time-reversed operator. Hence the 
expectation value of the pairing fields is (AJ,A„). This 



(AtA>=E 



2(2j + 1) 



(31) 



where rit = 



J im". 



are the neutron occupation num- 



bers. Any excess over the Fermi-gas value therefore in- 
dicates pairing correlations in the ground state. Our re- 
sults, which are plotted in Fig. 0] confirm a suppression 
of these correlations, as the contribution of the added 
neutrons to the pairing gradually decreases and the cor- 
relations become noticeably quenched beyond 96 Zr and 
100 Mo. 

In addition, the occupation numbers of various orbitals 
are plotted in Fig. demonstrating that additional neu- 
trons are distributed into the available orbitals rather 
uniformly while protons tend to migrate from the 0gg/2 
to the lpi/2 orbital. In the case of Mo isotopes, the 
spurious effect of exceeding the maximum- allowed occu- 
pancy for the lpi/2 orbital is a result of the extrapolation 
scheme that was used, and is indicative that the relative 
error bars on the occupation data after extrapolation are 
approximately 15% of the value of the occupation num- 
ber. Note that the deformation driving 0hn/ 2 remains 
only slightly occuppied. 
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FIG. 5: Orbital occupation numbers. 



IV. CONCLUSION AND PERSPECTIVES 

We have introduced a new approach for the implemen- 
tation of the Shell Model Monte Carlo method to perform 
shell-model calculations using non-identical proton and 
neutron valence spaces. General features of the SMMC 
method were reviewed. Differences between the isospin 
and the pn- formalisms have been pointed out; in partic- 
ular, the T z projection has been described in detail. 

The results of the SMMCpn approach have been val- 
idated in the sci-shell using a "good" -signed schematic 
interaction and in the /p-shell using the realistic KB3 
interaction. In the latter case, we dealt with the sign 
problem using an extrapolation method. 

As the first novel application of the SMMCpn ap- 
proach, we performed a set of calculations for the even- 
even 9 °- 104 Zr and 92_106 Mo isotopes, using a realistic ef- 
fective interaction in the valence space described by (it : 



lpi/2,059/2) and (v. ld 5 /2,2si/ 2 , ld3/2, 00 7 /2, Oftn/2) or- 
bitals. A comparison of the ground-state energies of the 
first few nuclei in both isotope chains showed excellent 
agreement with the exact diagonalization results and pro- 
vides a definitive test of our algorithm and the SMMCpn 
method. We then studied the transitional nature of the 
isotopes by using the B{E2) strength as a gross measure 
of the 0+ — 2+ separation. Along both isotope chains, we 
have obtained an enhancement in the B{E2) strengths 
as a function of the added neutrons, accompanied by a 
quenching in the neutron-pairing correlations. In spite 
of this qualitative reproduction of the onset of deforma- 
tions, further research clearly is needed for a qualitative 
result in the heavier Zr and Mo nuclei. A comparison 
with the experimental data suggests that this situation 
may be a shortcoming due to the degrees of freedom that 
are absent in the chosen valence space and possibly due 
to the value of the 0/i n / 2 orbital. We will investigate this 
further in future work. 

Apart from future applications involving realistic effec- 
tive interactions, use of schematic interactions in SMM- 
Cpn applications should be an interesting direction of 
research. Such interactions have been commonly used to 
calculate realistic estimates of collective properties and 
level densities; the latter is an important ingredient in 
the prediction of nuclear reaction rates in astrophysics. 
Parity dependence of these densities may play a crucial 
role in the nucleosynthesis. We believe that SMMCpn 
will prove to be a useful computational tool in this re- 
gard. 
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